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I. INTRODUCTION 

The front form [lj is widely recognized as an alterna- 
tive form of a relativistic theory. It has been successfully 
applied in particular in the context of quantum chromo 
dynamics (QCD), e.g., to describe deep inelastic scatter- 
ing or to model the structure of hadrons, see e.g. Ref. Q. 
The front form provides a rigorous relativistic many-body 
theory that has been applied, e.g., to investigate nuclear 
matter 0, 0] . Therefore it is quite natural to ask for the 
merits of light-front quantization in the context of quan- 
tum statistics to describe many-particle systems of finite 
temperature T and finite chemical potential fi. There is 
a need to do so, in particular, to model the evolution of 
the early universe, the dynamics of relativistic heavy ion 
collisions, multi-particle correlations in such a collision, 
relativistic Coulomb plasma, among others. 

The difficulty to generalize the statistical operator 
to the light-front is related to the longitudinal mode 
k + = k° + k 3 . Whereas temperature for the transverse 
modes related to k± — (k 1 , k 2 ) can easily be defined, the 
longitudinal modes, can only by related to the instant 
form by a limiting process (infinite momentum frame). 
Therefore the longitudinal modes appear in a proper def- 
inition of the statistical operator on the light-front. As a 
consequence the Fermi function (here given for a grand 
canonical ensemble) explicitly depends on the light-front 
energy k~ n , 



k 2 _+ m 2 
= fc+ 

and the light-front "mass" k + , viz. [f| 
f{k+,kj_)= expji (~Kn + \ k+ - l L 



(1) 



(2) 



This function is consistent with the instant form and has 
already been used to investigate the formation of three- 
quark clusters (nucleons) in the hot and dense phase of 
nuclear matter close to the QCD phase transition [f| . 

Recently, Brodsky 0, @ suggested a partition func- 
tion for a canonical ensemble (fj, = 0) that would lead 
to a different Fermi function (naive generalization). It 
has been explicitly shown that the naive generalization 
has serious problems, because the Fermi function leads 
to divergent integrals 9], which is not the case for the 
distribution function given in Eq. (J2J due to the presence 
of the k + term. Presently, there is a renewed interest and 
discussion on how to formulate thermal field theory on 
the light-front [13 EJ 13 

So far, none of these gen- 
eralizations consider finite chemical potentials that are 
relevant for QCD. Utilizing our previous results we gen- 
eralize the present formulation of thermal field theory on 
the light-front to include a finite chemical potential. 

As an example, we use the Nambu-Jona-Lasinio (NJL) 
model 0, that is a powerful tool to investigate the 
non-perturbative region of QCD as it exhibits sponta- 
neous breaking of chiral symmetry and the appearance 
of Goldstone bosons in a transparent way. The region of 
finite chemical potential and temperature can be investi- 
gated in relativistic heavy ion collisions of the newly con- 
structed (RHIC) and upcoming machines (LHC, GSI). It 
is particularly interesting since the phase structure is ex- 
pected to be complex ^(J El might be also 
relevant, e.g. for the physics of compact stars. Lattice 
QCD studies of this region (/j, ^ 0) are far more limited 
than for [i = [2(j, although some progress has been 
achieved recently, see e.g. [21| and refs. therein. 

In Sec.[n]we use our previous results to generalize the 
present formalism of thermal field theory on the light- 
front to finite chemical potentials. In Sec. IIIII we briefly 
introduce the NJL model on the light-front to give an 
example. This will be done along the lines of Ref. [2^ 
and in Sec. IIVI we will present numerical results for the 



2 



quark mass at finite temperature and densities (chiral 
restoration). 



II. LIGHT-FRONT STATISTICAL PHYSICS 

Here we present more details of our previous results, 
and provide the derivation of the distribution function 
for the grand canonical ensemble. See the Appendix [S] 
for notation. The four-momentum operator on the 
light-front is given by 



dS+ T +fl (x), 



(3) 



where T fJ,v (x) denotes the energy momentum tensor de- 
fined through the Lagrangian of the system and is the 
quantization surface, see Eq. (|A7|) . Hence the Hamilto- 
nian is given by 



dS+ T + -(x) 



(4) 



To investigate a grand canonical ensemble it is necessary 
to define the number operator. The number operator on 
the light-front is given by 



N 



J dS+ J + (x), 



(5) 



where J v {x) is the conserved current. These are the nec- 
essary ingredients to generalize the covariant calculations 
at finite temperature [UiHH to the light-front. The 
grand canonical partition operator on the light-front is 
given by 

Z G = cxp ( / dS+ [-p v T +v (x) + aJ+(x)]\ , (6) 



where a — fx/T, with the chemical potential p. The 
velocity of the medium is 



(7) 



with the normalized time-like vector u v u v — 1 [23|. We 
choose 



u v = (ti-,u+ tf- 1 ) = (1,1,0,0), 



(8) 



hence u + — u = 2u 3 = 0, i.e. the medium is at rest. 
The grand partition operator becomes 

Z G = e~ K/T , K = i(i>- + P+) - pN (9) 

with P ± and N defined in Eqs. and JSJ). To further 
evaluate Eq. © , the operators will be used in Fock space 



representation. For an ideal gas of Fermions we use the 
solution of the free Dirac field (conventions of [2|| ) 



dk+d 2 k ± 



b x (k)u aX (k)e-* k °- x 
+d{(k)v aX (k)e^ x ) (10) 



with the four- vector fc on = (fc on , k), the abbreviation k = 
(fc + , k±), and the Fock operators 

{b x (k),bUp)} = S(k+-p + )S^(k ± -p ± )5 X(T) (11) 
{d x (k),4(p)} = S{k + -p+)S^(k ± -p ± )S Xa . (12) 

The Hamiltonian takes the form 
P- = X) / rffc+rf2fc ^ fc on (b{(k)b x (k) + 4(S)d A (fc)) 

A ^ 

and for the plus component of the Dirac current, 
N = J2 [ dk+d 2 k ± (b[(k)b x (k) - d{(k)d x (k)y (14) 

Inserting all the above into © gives 

K = J2 I dk+d 2 k ± (V b{(k)b x (k) 

A ^ 

+K- d{(k)d X (k)), (15) 



k± = -{k- n + k+)T»- 



(16) 



The density operator for a grand canonical ensemble |2£ 
l27j in equilibrium follows 



PG = (Tr e-"/ J )-> e 



-K/Ts-l -K/T 



(17) 



The light-cone time-ordered Green function for 
fermions is 

G a p(x-y) = 0(x + ~y + ) g>p(x-y)+8(y+-x + ) G< p (x-y), 

(18) 

where 



GZpte-v) = -i(*i3(y)* a (x)), (19) 



We note here that the light-cone time-ordered Green 
function differs from the Feynman propagator Sp in the 
front form by a contact term (which is not the case in 
the instant form). Evaluating Eq. (|19f) for the vacuum, 
i.e. (...) = (0| . . . |0) (isolated case), the Green function 
is 



(2^) 4 



,—ik-ix—y) 



S F {k) 



2k+ ) ~ 



r 



G(k~,k) 



(20) 



where the Feynman propagator is given by 

"fk + m 



m 2 + ie 



(21) 



Eq. I|20|) coincides with the light-front propagator given 
previously in Ref. 28). Using the independent compo- 
nents instead of ^(x) in the definition of the chrono- 
logical Green function Eqs. H18fl and 11911 leads to a dif- 
ferent propagator suggested, e.g., in |9ll29l|. 

To evaluate the ensemble average (...) = Tr(po ■ ■ ■ ) 
of Eqs. I|18|) and i|19|) (in-medium case), we utilize the 
imaginary time formalism [26l [27j ■ We rotate the light- 



front time of the Green function to imaginary values, i.e. 
u • x — ► it. The anti-periodic boundary condition of the 
imaginary time Green function holds on the light-front 
as well. Hence the /c _ -integral is replaced by a sum of 
light- front Matsubara frequencies io n according to |5| , 



1 



1 



2 n ' 



(22) 



where w„ = 7rAT, A = 2n + 1 for fermions [A = 2n for 
bosons]. For noninteracting Dirac fields the imaginary 
time Green function becomes 



G{k-,k) 



jk 



8{V 



2 Kn 



2 "-on 



7^o 



-He 2k 
m 6(—k 



(1 -/+(*)) + T 



it" 

9 "<7l 



IE 



2k+ 



jk 



-f-(-k) + - 



kn 

jk 



Ik- 

2 "'on 



is 2k+ 
6{-k 



(23) 



2 ^on 



IE 



2k+ 



-(!-/-(-*))■ 



For a grand canonical ensemble the Fermi distribution 
functions of particles / + = / and antiparticles /~ are 
given by 



f ± (k + ,k 1 _) = 



e K± / T + 1 



<'xs><{ - Ufc on + -fc+TM I > + ■ 



(24) 



and fc~ n of Eq. The (particle) fermionic distribu- 

tion function on the light-front has been given in Eq. J3J). 
The propagator for this case (/~ = 0) has been given 
previously and used to investigate the stability and dis- 
sociation of a relativistic three-quark system in hot and 
dense quark matter (jj |(| . Note that for the canonical 
ensemble, fi = 0, this result coincides with the one given 
more recently in Q (up to different conventions in the 
metric) . 

For equilibrium the imaginary time formalism and the 
real time formalism are linked by the spectral func- 
tion mill. 



III. NJL MODEL ON THE LIGHT-FRONT 

The Nambu-Jona-Lasinio (NJL) originally suggested 
in 0, ^| has been reviewed in Ref. [3lJ as a model of 
quantum chromo dynamics (QCD), where also a general- 
ization to finite temperature and finite chemical potential 
has been discussed. Most of the applications of the NJL 
model have utilized the instant form quantization. On 
the other hand light-front quantization of quantum field 
theories has emerged as a promising method for solv- 
ing problems in the strong coupling regime. Examples 
and in particlar the light-front formulation of QCD, are 
reviewed in 0. Based on a Hamiltonian the light-front 
approach is particularly suited to treat bound states (cor- 
relations), the most relevant consisting of qq and qqq va- 
lence quarks. The formation of bound states is in partic- 
ular interesting in the vicinity of the quark-hadron phase 
transition. The front form is useful in the context of 
quantum statistics, since the Fock space representation 
allows a consistent formulation of thermodynamic prop- 
erties, as demonstrated in the previous Section. Since 
the NJL model, on one hand, is a powerful tool to in- 
vestigate the non-perturbative region of QCD and, on 
the other hand, is comparatively transparent and sim- 



4 



pie, we use it here as an example to tackle the following 
questions How are hadrons formed/dissociated? Is 
the formation/dissociation transition at the same place 
in the phase diagram as the chiral phase transition? Ap- 
proaches with a closer connection to QCD than the NJL 
model, based, e.g., on the method outlined in Ref. Q, 
have to follow. These may be of use as either complemen- 
tary or supplementary approaches to lattice QCD (e.g. at 
large chemical potentials /i, where lattice QCD is still not 
applicable). 

The NJL Lagrangian is given by 

C = 4>{ijd - m )^ + G ((^V) 2 + (V^75tV>) 2 ) ■ (25) 
In mean field approximation the gap equation is 

m = m -2G(ipip) = m + 2iGX [ -^-^ TrS F (k), (26) 

J (27T) 4 

where A = NfN c in Hartrcc and A = NfN c + ^ in 
Hartree-Fock approximation, N c (Nf) is the number of 
colors (flavors). 

A. Isolated case 

On the light-front the realization of chiral symmetry 
breaking is a subtle issue (3^, IH, 0- Unlike the 
instant form, for a given momentum (fc + ,fcjj, k + =/= 0, 
the dispersion relation on the light-front Q is unambigu- 
ously determined. This is usually denoted as the simple 
vacuum structure on the light-front. A subtlety is related 
to the zero modes, k + = 0, which contain the informa- 
tion on symmetry breaking. For a recent review on the 
realization of chiral symmetry breaking on the light-front 
see, e.g., Ref. H3- Here we basically follow Ref. p^ . 

After performing the fc~ integration the gap equation 
becomes (x > 0) 

- f dxd 2 k\ , x 

m-m^2GxJ^- f ^ (27) 

where in addition mo — > mo and G — > G have been renor- 
malized. The parameters mo, G and mo, G are related via 



The quark mass changes with the parameters T and fi. 
Note that being a single particle distribution the Fermi 
functions that appear in the gap equation depend on 
k + , rather than x = k + /P + . Therefore the LB regular- 
ization scheme cannot be used without modification that 



1/N C expansion (see appendix B of H^)- The explicit re- 
lation is given below for the regularization schemes dis- 
cussed here. Note that due to the trace either G{k) or 
S F (k) of Eq. (2} might be used, since Tr 7 + = 0. The 
integral of Eq. I|27|l clearly diverges and one has to in- 
troduce a regularization scheme. Several regularization 
schemes have been in use. A common scheme has been 
suggested by Lepage and Brodsky |3g (LB) 

k\ < Al B x{l -x)- m 2 , (28) 

which implies x± < x < x%, where 

x h2 = ~ (l T y/l 4mVA 2 B ^) . (29) 

The equivalence of the light front LB scheme to the in- 
stant form three momentum (3M) scheme, where (|26f) 
is integrated over k° and k 2 < A| M has been shown in 
Ref. [22 ] utilizing the Sawicki transformation |37j| . To 
compare results note that for the mentioned regulariza- 
tion schemes the following identifications hold 

m ,LB = "Z0,3M, (30) 
Glb = G3M, (31) 
Alb - 4(A 2 M +m 2 ). (32) 

Numerical values are given in the next section. The cal- 
culation of the pion mass m^, the pion decay constant 
/tt, and the condensate value follows (22| . 



B. In- medium case 

Here we are concerned in an application of the formal- 
ism to chiral restoration in hot and dense matter. The 
statistical approach shown in the previous section for fi- 
nite temperature T and chemical potential fi will be used. 
Because of the medium, the propagator to be used in the 
gap equation it26*|) is given in Eq. (|2*3*|l . The gap equation 
becomes 



(33) 

I 

will be shown below. To regularize Eq. I|33|l we require 
instead 



Kn + k+ < 2fi - (34) 



I 

/flk + fl 2 k\ 
2fc+ (27r )3 4m ( T ' M l - f + ( k+ > k ^ - f~( k+ > 
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As a consequence < k + < k^ and 

k\ < 2ttk + - {k+ f - to 2 , (35) 
fc+ 2 = ilT V 7 ^ 2 -m 2 . (36) 

To connect this fl scheme with the LB and the 3M 
schemes we use the Sawicki transformation [37I ]. 

k+ 1 / k 3 \ 

* = w = i{^h (37) 

UJ k = + \/fc 2 +(fc 3 ) 2 +TO 2 , (38) 

I 



where P + is evaluated in the rest system of the two par- 
ticles involved in the gap equation because of the loop, 
P + = 2ujk- Starting from Eq. (|33|l with the regulariza- 
tion conditions H35fl and (|36[) we substitute k + by x and 
replace the regularization conditions accordingly. If we 
choose 

A LB = 2f7, (39) 

we obtain the LB conditions l|28|l and Ij29(l . The resulting 
regularized gap equation in the LB scheme is 



m(T,fi) = to 0i lb + 2G L bA 



dxcPk± 



LB 



2x(2n 



3 4m(T,/i) (1 -f£ B (x,k±) -/ LB (a;,fcL)). 



f* B (x,k±) 



CX P i 



1 / k\ + TO 2 xA LB 



T \ 2xA 



LB 



(40) 

(41) 
(42) 



In the LB regularized gap equation the Fermi functions depend explictly on the regularization parameter Alb i because 
of the k + dependence of the Fermi function. Replacing k + by fc 3 in Eq. I|33|l we get 



m(T,fi) = m ,3M + 2G 3 mA 



d 3 k 



3M 



exp \ j, (uj k T M) r + 1 



2c; fe (27r)3 
n -1 



4m(T, M ) (l-/+(w fc ) -/"(a*)). 



(43) 
(44) 



and using (|35|l and (|36|l the regularisation condition fol- 
lows 

k 2 < fl 2 - to 2 = A 2 M . (45) 

This regularized version (|43() is analytically the same 
equation as given previously in Ref. [31| . Therefore we 
have shown explicitly that also in medium the SI regu- 
larized light front gap equation leads to the same results 
as the 3M regularized instant form gap equation. We re- 
mark here also that, because of the T and /i dependence 
of the mass, the regulators and Alb need to depend 
parametrically on temperature and chemical potential to 
keep the equivalence with the 3M scheme. 

IV. RESULTS 



mass too and the constituent mass is due to the finite 
condensate, which is (uu) 1 / 3 = —247 MeV for the pa- 
rameters given in Tabled The parameters are reasonably 
close to the cases used in the review Ref. |3l| . 

In hot and dense quark matter the surrounding 
medium leads to a change of the constituent quark mass 
due to the quasiparticle nature of the quark. The con- 
stituent mass as solution of Eq. (|33|l is plotted in Fig.^ as 
a function of temperature and chemical potential. The 
fall-off is related to chiral symmetry restoration, which 
would be complete for too = 0. It is related to the QCD 
phase transition. For T < 60 MeV the phase transition 
is first order, which is reflected by the steep change of the 
constituent mass. To keep close contact with the 3M re- 
sults we have chosen for the Q in-medium regulator mass 
fl 2 (T,n) = A 2 M + m 2 (T,n) with A 3M = 630 MeV fixed 



The model parameters are adjusted to the isolated 
system. We use the Hartree approximation, i.e. A = 
N c Nf — 6. Parameter values shown in Table |U are cho- 
sen to reproduce the pion mass m n = 140 MeV, the decay 
constant f n — 93 MeV, and to give a constituent quark 
mass of to = 336 MeV. The difference between the bare 



TABLE I: Parameters of the NJL model used in this analysis. 
Q, quoted for the isolated case. 



G[l(T 6 MeV] 


rh [MeV] 


n [MeV] 


5.51 


5.67 


714 



6 




o 



FIG. 1: Effective quark mass as a function of temperature 
and chemical potential. The fall-off is related to the vanishing 
condensate (uu), which shows the onset of chiral symmetry 
restoration. Critical temperature at fi = is T c ~ 190 MeV. 

250| J 




chemical potential ufMeV] 

FIG. 2: Chiral phase transition as defined in p5ftl |. The lower 
part is the chiral broken phase, whereas the upper part reflects 
the restored phase. 



suiting Fermi function depends on transverse and also 
on the k + momentum components. The k + compo- 
nents emerge in a natural way in a covariant approach 
and are also essential to fulfill the light front analog of 
the Thouless criterion [4C| for the appearance of Cooper 
poles (color superconductivity) As an application 

we have revisited the NJL low energy model of QCD. 
Because of the zero range interaction a regularization is 
necessary. The equivalence of instant form 3M regular- 
ization and light front Q regularization utilized here for 
the gap equation holds for a T and /i dependent cut-off 
Q 2 (T,n) = A| M + m 2 (T 7 fi). Hence we reproduce the 
phenomenology of the the NJL model, in particular the 
gap-equation and the chiral phase transition. The NJL 
model constitutes a nontrivial example of the equivalence 
between instant form and front form quantization. We 
argue that the light front formulation of thermal field 
theory is a very useful and powerful alternative tool to 
make closer connection to QCD, in particular for the light 
front formulation of QCD, e.g., along the lines of The 
region of interest is that of finite //, where lattice QCD 
just begins to become available for small /i 21]. Fur- 
ther merits of the light front approach are obvious when 
considering multi-quark correlations as has been already 
outlined in 
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APPENDIX A: LIGHT-FRONT NOTATION 



for all T and fi. 

We define the phase transition to occur at a tempera- 
ture at which m(T, fi) is half of the isolated constituent 
quark mass [3c| . Because of the above mentioned equiv- 
alence we do not expect any qualitative difference to the 
results of the instant form analysis provided, e.g. by 
Refs. [3l|, |3!|. The phase diagram is shown in Fig. [21 
The line indicates the phase transition between the chi- 
ral broken and the unbroken phase. 

V. CONCLUSION 

We have given a relativistic formulation of thermal field 
theory utilizing the light front form. The proper parti- 
tion operator (and the statistical operator) have been 
given for the grand canonical ensemble. The special case 
of a canonical ensemble is given for \i = 0. The re- 



The contravariant components of a four vector in light- 
front coordinates may be written as 

a M = (a~ , a + : a 1 , a 2 ) = (a~,a + ,a^) = (a - , a) (Al) 

If we define a ± = a ± a 3 , the standard Lorentz scalar 
product, can be recovered by defining the following met- 
ric tensor 



9^ = 



(0 


1 

2 





o ^ 






2 





o \ 


1 

2 











, <r = 


2 

















-1 











-1 





\0 












{0 








-1/ 



(A2) 
i.e. 

OuV = g^b v = \{a-b+ + o+6-) - a x 6±. (A3) 
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Hence the covariant components of the four vector are 
given by 

a M = (a-,a + ,ai,a 2 ) = (§a + , \cC , -ax) = (a_,a). 

(A4) 

Note 

1 , 1 /4 . 

a- = 2 a > a + = 2 a ' ( 5 ) 

The four-dimensional volume element is 

d 4 x = dx°dx 1 dx 2 dx 3 = ^dx~dx + d 2 x± (A6) 

The three-dimensional volume element on the light-like 
surface S v characterized by the light-like vector ui v in the 
plane with lli v x v = x + = const, is 

d i x\ a = —dx~ d 2 x± = dS+. (A7) 
16+2 "' 



The scalar product of the light-like surface element dS, 
with a vector a v is then 



dS v a v = -dx'd 2 x ± a+. (A8) 



From the 7-matrices 7 M = (7 + ,7 , 7x) it is possible to 
define Hermitian projection operators 



A ± = i = 7 T 7 ± - (A9) 
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